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ABSTRACT 

Gravitational stability of a disc consisting of the gaseous and the stellar components 
are studied in the linear regime when the gaseous component is turbulent. A phe- 
nomenological approach is adopted to describe the turbulence, in which not only the 
effective surface density but the velocity dispersion of the gaseous component are both 
scale-dependent as power-law functions of the wavenumber of the perturbations. Also, 
the stellar component which has gravitational interaction with the gas is considered as 
a fluid. We calculate growth rate of the perturbations and in the most of the studied 
cases, the stability of the disc highly depends on the existence of the stars and the 
exponents of the functions for describing the turbulence. Our analysis suggests that 
the conventional gas and star threshold is not adequate for analyzing stability of the 
two-component discs when turbulence is considered. 
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1 INTRODUCTION 

Various instability mechanisms may operate in an accretion 
disc depending on the physical properties of the disc itself. 
An imbalance of the heating and cooling rates in an thin ac- 
creti on disc may lead to a disc which is thermally unstabl e 
(e.g., lLightman fc Eardlevlll974l ; |Pringlej|l9761 : |Piranlll978h . 
Under some conditions, an accretion disc is hot enough to 
be fully or partially ionized which then implies to expect sig- 
nificant dynamical effects of the magnetic field. In fact, an 
instability related to the rotational profile and the magnetic 
field of an accretion disc is known to be the driven mech- 
anism of generating tur bulence inside the accretion discs 
|Balbus fc Hawlevlll99ll ). In the most of the accreting sys- 
tems, the mass of the disc is negligible in comparison to the 
mass of the central object and so, one may safely neglect 
the self-gravity of the disc. However, gaseous or stellar discs 
are subject to the gravitatio nal instability if the Toomre pa- 
rameter drops below unity (|Toomrelll964l ') and the cooling 
time scale is less th an a few time of the dynamical time 
scale (G ammic 2001). Form ation of structures like planet s 
in protoplanetary discs (e.g., Boss 1998; Rafikov 2005. 



or stars at the Galactic center (e.g., Goodman & Tan 



2007) 



iNavakshin fe Cuadrdl20^ : lNavakshinll2006l : lLevinll2007l ~ 

described based on the gravitational instability. 



2004 
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Analysis of the gr avitationa l stabi lity of a disc starte d 
by the pioneer works of Safronov ( 1960) and Toomre ( 1964). 
The onset of axisymmetric gravita tional instabil ity is deter- 
mined by the following condition (|Toomrelll964l ) 



Q = ^k<Qo, (i) 

where Qo = 1 and c s , k and E are the sound speed, epicyclic 
frequency and the surface density of the disc, respectively. 
The non-dimensional Toomre parameter Q has a vital role 
in the gravitational stability of an accretion disc. The non- 
axisymmetric perturbations grow for somewhat higher Qo- 
However, the Toomre condition (fTJ is not the only 
criteria for the fragmentation of a disc. In fact, thermo- 
dynamics of the disc and its efficiency of cooling dictates 
another condition for the fragmentat ion of a disc a s has 
been suggested in a seminal work by iGammi c (2001) and 
investigated via ex te nsive nu merical simulations (e.g., 



2005: 



Johnson fc Gammid 120031: iRice. Lodato fc Arrnitagel 



Clarke. Harper-Clark fc Lodatol 



Cossins. Lodato fc Clarkdl2009l) 



2007; 



Although validity of the Toomre condition for the sta- 
bility of the self-gravitating discs has also been confirmed 
by direct numerical simulations, in some of the astrophys- 
ical systems like galaxies one component description is not 
adequate to study t h eir stability (e.g ., lYim et af I 1201 ll : 
iBournaud et all l2010l : lElmegreenl |20 1 ll ) . A typical galaxy 
consists of the gas and the stars. Over recent two decades, 
many authors studied gravitational stability of such two- 
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components discs (e.g. , Jog fc Solomon! Il 9841 ; lRomeolll992l ; 
IWang fc SiUdll994IJoelll996l ). The main outcome is the stel- 
lar component has a destabilizing role either a fluid approx- 
imation or a collisionless description is used to model the 
stellar component. Unfortunately, it is not possible to rep- 
resent the stability condition of a two-component disc is a 
closed analytical form like the one component case, though 
some authors introduced approximate analyt ical formulae 
for the onset of the instability in a disc (e.g. JWang fc Silkl 
ll994lRomeo fc Wiegertll201lf) 



Recently. iRomeo. Burkert fc Agert j (|20T0h (hereafter; 
RBA) studied possible effects of the turbulence on the 
gravitational stability of the accretion discs using a sim- 
ple, but very illustrative model for the turbulence. They 
used scale-dependent relations for the effective surface den- 
sity and the velocity dispersion of the gaseous component. 
Then, they could obtain conditions under which the self- 
gravitating disc becomes unstable depending on the prop- 
erties of the turbulence. For a range of the input param- 
eters, they showed the conventional Toomre parameter is 
modified to correctly describe the stability of the disc when 
turbulence of the gas is considered. However, such kind of 
approach for studying gravitational stability of the discs 
has already been applied by a few authors, but for some 
specifi c cases (e.g., Elm egreen 1996). In another relevant 
study. lElmegreenT(|201ll ) investigated gravitational instabil- 
ities in two-component disc galaxies when the gas dissipates 
on the local crossing time. However, he did the linear anal- 
ysis by considering the energy equation and prescribing a 
scale-dependent dissipation rate instead of following a phe- 
nomenogical approach to describe the turbulence. The main 
finding is an increase to the conventional gas and star thresh- 
old by a factor of two or three because of the existence of the 
dissipation which removes the stabilization from turbulent 
pressure. 

In this paper, we study gravitational stability of a two 
fluid disc where the turbulence of the gaseous component 
is considered similar to RBA. However, we will show that 
different regimes of the (in)stability according to RBA are 
significantly modified in t he presence of the star s. One can 
also compare our results to lJog fc Solomon! l|l984h to see how 
stability of a two fluid disc is changed when the turbulence 
of the gaseous component is considered. In the next section, 
our basic assumptions and the main dispersion relation are 
presented. Then, we will analyze stability of the system for 
a wide range of the input parameters and the results are 
compared to RBA. 



2 LINEAR PERTURBATIONS 

The basi c equations of ou r mode l are identical to that intro- 
duced in IJog fc Solomon! l|l984h . in particular regarding to 
the two-fluid description of the system. Gaseous and the stel- 
lar components are considered as two separate fluids which 
can only interact gravitationally. Although we assume the 
stellar component is collisional (i.e., fluid description), it is 
also possible to consider the c ollisionless nature of the stel- 
lar component |Rafikovll200lh . Both the gaseous and the 
stellar components are assumed to be isothermal character- 
ized by the sound speed a s and the velocity dispersion er s , 



respectively. The unperturbed surface densities for the gas 
and the stellar fluid are denoted by E g o and E s o- Also, the 
scale-heights of the gaseous component and the stellar mass 
distribution are given by 2/i g and 2h s . But, we include turbu- 
lence of the gaseous component in a way that was introduced 
by RBA. In this approach, not only the initial density but 
also velocity dispersion of the gaseous component are both 
scale-dependent because of the existence of the turbulence: 



(2) 
(3) 



where a and b are input parameters which describe the na- 
ture of the turbulence. Density fluctuations has a power-law 
spe ctrum, E (k) oc k~ r , which implies p oc fc - ^ -1 ^ 2 (see, 
e.g. lElmegreen fc Scaloll2004 IScalo fc Elmegreenll2004 ). As- 
suming E e ff ~ ph g , one can easily show that a — (l/2)(r— 1). 
Also, the spectrum of the velocity fluctuations is power- 
law, i.e. E v (k) oc k~ s , and then a s oc fc _ ' s_1 ' / ' 2 . The above 
power-law relations are appropriate for the cold interstellar 
medium and the length 1/fco is the fiducial scale at which 
Toomre parameter and other stability quantities are mea- 
sured (RBA). Moreover, the mass-size scaling relation im- 
plies —2 a ^ 1 (RBA). However, observational and the- 
oretical studies imply more constraints on the acceptable 
ranges for a and b. HI observations show that a Kolmogorov 
spectrum can describe the density and the velocity fluctu- 
ations in these regions, i.e. a ~ 1/3 for the scales less than 
10 kpc, and b ~ 1/3 for the scales less tha n 1 kpc (e.g., 
lLazarian fc Pogosvan|[200ol ; iDutta et alj|2009t) . But in giant 
molecular clouds, the dens ity and the ve l ocity exponents ar e 
a ~ and b ~ 1/2 (e.g., iLarsonl Il98l"l ; iHever et alj|2009h . 
On the other hand, high and low resolution simulations of 
supersonic turbulence imply the typical values of (a, 6) as 
(1/2,1/2) and (2/3,1/2), respectively (RBA). 

If we assume a functional form as exp[i(kr + cot)] for the 
perturbed variables in the linearized hydrodynamical equa- 
tions, the dispersion relation is obtained as 



■to (a B + « g ) + (a B a g - /3 s /3 g ) = 0, 



where 



2 . , 2 2 

q s = k + k a s 



- 2nGkE s0 [1 - exp{-kh s )/kh s ] , 
Q g = k 2 + fc 2 f 2 — 2nGkT, c f{ [1 — exp(—khg)/kh g ] . 
P s = 27rGfcE s0 [1 - exp(-kh s )/kh s ] , 
/3 g = 27rGfcE cff [1 - exp(-kh g )/kh g ] . 



(4) 

(5) 
(6) 
(7) 
(8) 



We identify uj as the growth rate of the gas and star grav- 
itational instability. Since only the two-dimensional forcing 
is considered, a correction factor of [1 — exp ( — kh B )\/(kh s ) 
appe ared in the above equations (|Toomre! 1 1964 ). How- 
ever, IVandervoortl (|l970h introduced a correction factor of 
(1 + fc/ts)" 1 from a more detailed analysis. Note that the ap- 
propriate correction factor appears for the gas as well. Thus, 
if we set a = b = 0, it means that the correction factors are 
considered for the stellar fluid and the non-turbulent gas 
component and the equations are reduced to JS. Equations 
is solved analytically as 
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Figure 1. This Figure shows uj 2 versus A when a = 1 and b ^ 
1 and fco = 0.5fcx (top), fco = Ait (middle) and fco = l-5fcx 
(bottom). The other input parameters are k = 65 km s —1 kpc -1 , 
a s = 34.7 km s _1 , cr g0 = 5 km s _1 , E g o = 19 Mq pc -2 and 
Sso = 190 M© pc -2 . Dashed curve corresponds to a case without 
turbulence. 
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(a s + ag) ± y/(a s + Q g) 2 " 4 K Q g " frA) .(9) 



where only the following root leads to the instability (JS), 
1 



uT(fc) 



(a s + a s ) - yj (a B + a g ) 2 - 4(a s a g - /3 B /3 g ) 



(10) 



We note that if the stellar or the gaseous contributions to the 
(in)stability of the system are neglected, then growth rate of 
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the unstable perturbations is determined based on the sign 
of q s or Qf g . In other words, the stellar or the gaseous com- 
ponents are unstable when a s < or a g < 0, respectively. 

For a purely gaseous disc, RBA studied the instabil- 
ity regimes corresponding to a g < for different values of 
a and b. Now, we illustrate how this parameter study is 
modified when the stellar contribution is considered. Thus, 



we must determine under what circumstances equation (|10|l 
gives a negative value, i.e. uj 2 < 0. However, it is very un- 
likely to present the stability condition based on the growth 
rate equation (|10[1 in a closed analytical form. Instead, we 
adopt the parameters corresponding to the Galaxy as an il- 
lustrative example, but for different values of a and b. For 
a gaseous disc, RBA showed that the stability of the disc is 
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classified into seven categories depending on the values of a 
and 6. Now, we explore how these different regimes of the 
instability are modified in the presence of a stellar fluid. 



3 ANALYSIS 

JS adopted their input parameters for the Galaxy based on 
the observational data. We also consider a similar set of the 
input parameters, because it helps us to compare our results 
to the previous findings. In all plots, our input parameters 
correspond approximately to conditions at the distance R = 
6 kpc from center of the Galaxy. Although the stars in the 
disc have a range of velocity dispersions, we characterize 
the stellar fluid by a single velocity dispersion a s = 34.7 km 
s . The epicyclic frequency is k = 65 km s _1 kpc -1 . The 
sound speed in the gas and the ratio of the gas density to 
the stellar density are 5 km s _1 and 0.1, respectively. Also, 
we assume E g0 = 19 M pc -2 and E s0 = 190 M @ pc -2 . We 
also assume the disc is in vertical equilibrium and the scale 
heights are assumed 2h s » 150 pc and 2h s ~ 180 pc (JS). 

Having the above input parameters, we can plot growth 
rate versus the wavenumber of the perturbations using equa- 
tion (|10|l . Figure [T] shows up" versus the wavelength of the 
perturbations, A, when a = 1 and b ^ 1. In this case, RBA 
showed that a gaseous disc is stable for all wavenumbers as 
long as fco ^ &t = K 2 /(27tGEgo). Here, fcx is the Toomre 
wavenumber. But we found this condition is not valid when 
the stellar fluid is considered. In this case, we note that for 
the gaseous component the self-gravity term is independent 
of the wavenumber, but for the stellar fluid it is not. The 
Kolmogorov turbulence corresponds to a = 1 and 6=1/3 
which lies in the explored regime of Figure [1] HI observa- 
tions also suggest a Kolmogoro v scaling for the velocity and 
the density fluctu ations (e.g., iLazarian fc Pogosvar] |2000| ; 
iDutta et al.ll2009h . The growth rate curves in Figure [1] are 
drawn for three values of fco, i-e. fco = 0.5fcT (top plot), 
fco = fcT (middle plot) and fco = 1.5fcr (bottom plot). In 



all plots, dashed curves correspond to a case without tur- 
bulence, i.e. a = b = 0. Top plot of Figure [T] clearly shows 
the system is not stable for all the wavenumbers and always 
there is a range of fc where the disc is unstable, irrespec- 
tive of the value of b. But when the critical wavenumber 
fco increases, more interesting cases emerge according to the 
middle and bottom plots of Figure [1] For example, RBA 
predicts that for fco = 1.5fcT a gaseous disc is always unsta- 
ble, but bottom plot of Figure [T] shows a stability behavior 
when 6 = 2. Also, the range of fc for which the disc is unsta- 
ble becomes larger when fco increases. 

In Figure [3] we explore another case: a ^ 1 and 6 = 
0.5(1 + a). So, both the pressure and the self-gravity terms 
of the gaseous component have the same fc— dependence. 
Each curve is labeled by the corresponding values of a and 
6 and again the dashed curves show growth rate of the 
perturbations when the turbulence of the gas is neglected. 
Three different values of fco are adopted. For this case, RBA 
showed that a gaseous disc is stable for all fc as log as 
fco *5 fcj = (27rGEgo)/a- 2 , where fcj is the conventional Jeans 
wavenumber. However, we see that this condition is modi- 
fied at least for some of the input parameters, e.g. a = and 
6 = 0.5. 

Figures [3] and U show growth rate of the perturbations 
for the case a < 1 and 6 > 0.5(1 + a). RBA did not explore 
this case in detail because the zero(s) of w 2 (fc) should be de- 
termined numerically, case by case. But they mentioned one 
component gaseous disc is unstable at small scales, irrespec- 
tive of the existence of the turbulence. Here, we also study 
this case for a limited range of the input parameters. For a 
particular s et of t he parameters a and 6 from this category, 
Elmegrcen (|l996f ) has also studied stability of a turbulent 
gaseous disc corresponding to the Larson-type scaling rela- 
tions, i.e. a — — 1 and 6 = 1/2. He found that the disc is 
always stable at large scales and unstable at small scales. In 
both Figures [3] and 21 different values of fco are considered. 
We can see the disc becomes more stable with increasing fco 
so that when fco = 1.5fcj the system is stable for the whole 
values of a and 6. But when fco = 0.2fcj and a = 0, the sys- 
tem is unstable at two intervals of the wavenumber. These 
intervals of the instability diminish with increasing fco . Com- 
paring Figures [3] and [4] we can say the system becomes more 
stable as the value of a decreases. So, not only the results of 
RBA but Elmegreen's result are modified in the presence of 
a stellar fluid. 

Figure [5] shows growth rate of the perturbations for an- 
other considered case by RBA, i.e. a > 1 and 6 < 0.5(1 + a). 
For this case, RBA showed that a turbulent gaseous disc 
is unstable at large scale. Our Figure [5] confirms this be- 
haviour for a turbulent gaseous disc in the presence of the 
stellar fluid. 

The case a < 1 and 6 < 0.5(1 + a) is a Toomre-like sit- 
uation as discussed by RBA. They found a stability thresh- 
old Q in terms of the input parameters and the turbulent 
gaseous disc is stable at all wavenumbers if Q ^ Q. Figure 
[S] explores this case. Each curve is labeled by the relation 
between the corresponding parameters Q and Q based on 
RBA. For example, curves corresponding to fco = 1.5fcT and 
fco = 0.75/ct in the middle plot show stability in contrary to 
inequality Q < Q. In other plots, the system becomes un- 
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stable for the limited ranges of the wavenumber of the per- 
turbations. In all the plots, we also note that for fixed E s o 
and E g o, varying fco has the effect of varying the effective 
surface density. Thus, the system has higher effective sur- 
face density and becomes more unstable at lower fco when 
a < 0. But for positive values of a, we have less effective 
surface density and more stability at lower fco- 



4 CONCLUSIONS 

We studied gravitational stability of a disc consisting of two 
components (i.e., gaseous and stellar fluids) using a per- 
turbation analysis. We also considered the turbulent nature 
of the gaseous fluid using a phenomenological description 
which is supported by the observations and the numerical 
simulations. Although stab ility of a two-fluid disc without 
turbulence of the gas (e.g.Jjog fc Solomor]|l984j ; |jodll99^ ; 
IWang fc S ilk 1994; Rafi kovll200ir) or one compone nt turbu- 
lent gaseous disc ( Romeo. Burkert Sz Agert3l2010| ) have al- 
ready been studied analytically, but our analysis is a first 
step towards considering not only the two-fluid nature of 
such systems but the turbulence of the gas. We found that 
turbulence of gas highly affects the stability of the two- 
component stellar and fluid discs. According to our anal- 
ysis, not only there is not a closed analytical relation for 
onset of the instability in a two-component disc, but the 
conventional gas and star threshold does not correctly ad- 
dress the stability of such systems in the presence of tur- 
bulence of the gas. We also explored situations without the 
thickness correction for the gas. In fact, this correction fac- 
tor has a stabilizing effect, in particular at high wavenum- 
ber k. For the perturbations much larger than the thickness 
of the disc, the thickness correction becomes unity. But as 
the wavenumber increases so that for the scales larger than 
the thickness of the disc, the correction factor introduces an 
additional inverse k dependence to the effective surface den- 
sity. This weakens the self-gravitational effect significantly. 
In the presence of the turbulence of the gas, we also found 
this stabilizing effect. 

In dwarf galaxies or the far-outer regions of the galaxy 
discs, the conve ntional stability criteria suggests a high level 
of stabi lity (e.g.. Ivan Zee et al.lil997l ; lHunter et al.ll201ll ). Al- 



though [Elmegreen (2011) showed that such regions are more 
unstable when turbulence of the gas is considered, our anal- 
ysis shows that the stability of such systems highly depends 
on the nature of the turbulence (i.e., values of a and b) in 
the presence of stars. However, detailed numerical simula- 
tions of the galaxy discs (including the turbulence of the 
gas component) are needed to confirm our linear stability 
analysis. 
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